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Abstract. We discuss the leading particle spectra and diffractive mass 
spectra from the novel point of view, namely by treating them as partic- 
ular examples of the general energy flow phenomena taking place in the 
multiparticle production processes. We argue that they show a high de- 
gree of universality what allows for their simple description in terms of the 
Interacting Gluon Model developed by us some time ago. 



1. Introduction 

The multiparticle production processes are the most complicated phenom- 
ena as far as the number of finally involved degrees of freedom is concerned. 
They are also comprising the bulk of all inelastic collisions and therefore are 
very important - if not per se than as a possible background to some other, 
more specialized reactions measured at high energy collisions of different 
kinds. The high number of degrees of freedom calls inevitably for some kind 
of statistical descrition when addressing such processes. However, all corre- 
sponding models have to be supplemented by information on the fraction 
of the initial energy deposited in the initial object (s) (like fireball(s)) being 
then the subject of further investigations. 

Some time ago we have developed a model describing such energy de- 
posit (known sometimes as inelasticity) and connecting it with the apparent 
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dominance of multiparticle production processes by the gluonic content of 
the impinging hadrons, hence its name: Interacting Gluon Model (IGM) 
(Fowler et al, 1989). Its classical application to description of inelasticity 
(Duraes et al, 1993) and multiparticle production processes in hydrody- 
namical model approach (Duraes et al, 1994) was soon followed by more 
refined applications to the leading charm production (Duraes et al, 1995) 
and to the (single) diffraction dissociation, both in hadronic reactions (Du- 
raes et al, 1997a) and in reactions originated by photons (Duraes et al., 
1997b). These works allowed for providing the systematic description of 
the leading particle spectra (which turned out to be very sensitive to the 
presence of diffractive component in the calculations, not accounted for be- 
fore) (Duraes et al, 1998a) and clearly demonstrated that they are very 
sensitive to the amount of gluonic component in the diffracted hadron (Du- 
raes et al, 1998b) and (Duraes et al, 1998a). We have found it amusing 
that all the results above were obtained using the same set of basic parame- 
ters with differences arising essentially only because of different kinematical 
limits present in each particular application (i.e., in different allowed phase 
space). All this points towards the kind of universality of energy flow pat- 
terns in all the above mentioned reactions. 

Two recent developments prompted us to return again to the IGM ideas 
of energy flow: one was connected with the new, more refined data on the 
leading proton spectra in ep — > e'pX obtained recently by ZEUS collab- 
oration ^ (which are apparently different from what has been used by us 
before in (Diiracs et al., 1998a), (Duraes et al, 1998b) and (Duraes et 
al, 1998a)). The other was recent work on the central mass production in 
Double Pomeron Exchange (DPE) process reported in (Brandt et al., 2002) 
allowing in principle for deduction of the Pomeron-Pomeron total cross sec- 
tion a IP- IP . In what follows we shall therefore provide a brief description 
of IGM, stressing the universality of energy flow it provides and illustrating 
it by some selected examples from our previous works. The new results of 
ZEUS will be then shown again and commented. Finally, we shall present 
our recent application of the IGM to the DPE processes as well (Duraes et 
al, 2002). 

2. IGM and some of its earlier applications 

The main idea of the model is that nucleon-nucleon collisions (or any 
hadronic collisions in general) at high energies can be treated as an in- 
coherent sum of multiple gluon-gluon collisions, the valence quarks playing 
a secondary role in particle production. While this idea is well accepted for 

^Private information from A.Garfagnini, see also ZEUS Collab. presentation in these 
proceedings. 
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large momentum transfer between the colliding partons, being on the basis 
of some models of minijet and jet production (for example HIJING (Wang 
et al, 92)), in the IGM its validity is extended down to low momentum 
transfers, only slightly larger than A.qcd- At first sight this is not justi- 
fied because at lower scales there are no independent gluons, but rather 
a highly correlated configuration of color fields. There are, however, some 
indications coming from lattice QCD calculations, that these soft gluon 
modes are not so strongly correlated. One of them is the result obtained 
in (Giacomo et al, 1992), namely that the typical correlation length of the 
soft gluon fields is close to 0.3 fm. Since this length is still much smaller 
than the typical hadron size, the gluon fields can, in a first approximation, 
be treated as uncorrelated. Another independent result concerns the de- 
termination of the typical instanton size in the QCD vaccum, which turns 
out to be of the order of 0.3 fm (Shacfcr et al, 1998). As it is well known 
(and has been recently applied to high energy nucleon-nucleon and nucleus- 
nucleus collisions) instantons are very important as mediators of soft gluon 
interactions (Shacfcr, 1998). The small size of the active instantons lead 
to short distance interactions between soft gluons, which can be treated as 
independent. 

These two results taken together suggest that a collision between the 
two gluon clouds (surrounding the valence quarks) may be viewed as a sum 
of independent binary gluon-gluon collisions, which is the basic idea of our 
model. The interaction follows then the usual IGM picture (Fowler et ai, 
1989) and (Duraes et al, 1998a), namely: the valence quarks fly through 
essentially undisturbed whereas the gluonic clouds of both projectiles in- 
teract strongly with each other (by gluonic clouds we understand a sort of 
"effective gluons" which include also their fluctuations seen as qq sea pairs) 
forming a kind of central fireball (CF) of mass M. The two impinging pro- 
jectiles (usually protons/antiprotons and mesons) loose fractions x and y 
of their original momenta and get excited forming what we call leading jets 
(LJ's) carrying Xp = 1 — x and Xp = 1 — y fractions of the initial momenta. 
Depending on the type of the process under consideration one encounters 
different situation depicted in Fig. 1. In non-diffractivc (ND) processes one 
is mainly interested only in CF of mass M, in single diffractive (SD) ones 
in masses Mx or My (comprising also the mass of CF) whereas in double 
Pomeron exchanges (DPE) in a special kind of CF of mass Mxy- The only 
difference between ND and SD or DPE processes is that in the later ones 
the energy deposition is done by a restricted bunch of gluons which in our 
language are forming what is regarded as a kind of "kinematical" Pomeron 
(IP), the name which we shall use in what follows. 

The central quantity in IGM is then the probability to form a CF car- 
rying momentum fractions x and y of two colliding hadrons (Fowler et al., 
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Figure 1. Left panel: schematic IGM pictures for (a) non-diffractive (ND), (&) and 
(c) single diffractive (SD) and (d) double Pomeron exchange (DPE) processes. Their 
corresponding phase space limits are displayed on the right panel. The | In - line in the 
right (a) panel indicates the rapidity Y of the produced mass M . 



1989) and (Duraes et al, 1998a) which is given by: 
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where i = S,H, and = Kmin whereas K!^^^ = 2;>t„.„/i/s. The values 

of Xmax and Umax depend on the type of the process under consideration 
(cf. Fig. 2). For ND processes Xmax = Umax = 1 (all phase space above the 
minimal one is allowed) whereas for SD and DPE there are limitations seen 
in Fig. 2 and discussed in more detail in the appropriate sections below. 
Here G's denote the effective number of gluons from the corresponding pro- 
jectiles (approximated in all our works by the respective gluonic structure 
functions), and are the soft and semihard gluonic cross sections, 
PTmin minimum transverse momentum for minijet production and a 

denotes the impinging projectiles cross section. 

Lets us close this section by mentioning that, as has been shown in 
(Fowler et al, 1989), (Duraes et al, 1994) and (Duraes et al, 1995), IGM 
can describe both the hadronic and nuclear collision data (providing initial 
conditions for the Landau Hydrodynamical Model of hadronization (Duraes 
et al, 1994)) as well as some peculiar, apparently unexpected features in the 
leading charm production (Duraes et al, 1995) (mainly via its strict energy- 
conservation introducing strong correlations between production from the 
CF and LJ). It was done with the same form of the gluonic structure 
function in the nucleon used: G{x) = p{m + 1)(1 — x)'"'' /x with defold 
value of m = 5 and with the fraction of the energy-momentum allocated to 
gluons equal to p = 0.5 and with a^gg {xys) = const (notice that results are 
sensitive only to the combination of p^Ggg/a). 

3. Single DifFractive processes in the IGM 

In the last years, diffractice scattering processes received increasing atten- 
tion mainly because of their potential ability to provide information about 
the most important object in the Regge theory, namely the Pomeron (IP), 
its quark-gluonic structure and cross sections. Not entering the whole dis- 
cussion (Goulianos, 1983) ^ we would like to show here the possible ap- 
proach, based on the IGM, towards the mass(Mx) distributions provided 
by different experiments. In Figs, lb and Ic the understanding of what 
such mass means in terms if the IGM is clearly shown: it contains both the 
central fireball and the LJ from the initial particle which got excited. The 
only difference between it and, say, the corresponding object which could 
be formed also in Fig. la is that the energy transfer from the diffracted 
projectile is now done by the highly correlated bunch of gluons dcnnoted 
IP which are supposed to be in the colour singlet state. The other fea- 
ture also seen in Fig. 1 (right panel) is that now only the limited part of 



^See also talk by K. Goulianos in these proceedings. 
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the phase space supporting the xi^^ u) distribution is ahowed and that the 
Umits depend on the mass Mx we are going to produce (and observe). 
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Figure 2. Examples of diffractive mass spectra for (a) pp collisions compared with 
Tevatron data (Abe et al, 1994) (Fig. 4 from (Duraes et al, 1997a)) and for (b) -yp 
collisions compared with HI data (Adloff et al, 1997) (Fig. 2b from (Duraes et al, 
1997b), Vector Dominance Model has been used here with G''{x) = (x) = 2p{l — x)/x; 
the different curves correspond to different choices of (mo GeV, a mb ): 1= (0.31,2.7), 
11= (0.35,2.7), 111= (0.31,5.4) and IV= (0.35,5.4)). 



In technical terms it means that in comparison to the previous appU- 
cations of the IGM we are free to change both the possible shape of the 
function Gip{x) (telling us the number of gluons participating in the pro- 
cess) and the cross section a in the spectral function uj in eq. (4) above. 
Actually we have found that we can keep the shape of G{x) the same as 
before and the only change necessary (and sufficient) to reach agreement 
with data is the amount of energy-momentum p = pjp allocated to the im- 
pinging hadron and which will find its way in the object we call IP. It turns 
out that Pip ~ 0.05 (to be compared with p ~ 0.5 for all gluons encountered 
so far). In Fig. 3 we provide a sample of results taken from (Duraes et al., 
1997a) and (Duraes et al, 1997b). They all have been obtained by putting 
Xmax = 1 and Umax = / s in eq. (2) above and by writing 

dM^ = <^y {^x - sy) [xy - Kl,i^ . (5) 

As can be explicitely shown the characteristic behaviour of diffractive 

mass spectra are due to the G{x) ~ 1/x behaviour of the gluonic structure 
functions for small x. The full formula results in small deviations following 
precisely the trend provided by experimental data (and usually attributed 
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in the Regge model approach to the presence of additional Reggeons (Gou- 
lianos, 1983)). 

4. Leading Particle spectra in the IGM 

With the above development of the IGM one can now think about the 
systematic survey of the leading particle spectra, both in hadronic and in 
7P collisions (Duraes et al, 1998a) (cf. Figs. 4 and 5). The specific prediction 
of the IGM connected with the amount of gluons in the hadron available 
for interactions (i.e., for the slowing down of the original quark content of 
the projectile) has been discussed in (Duraes et ai, 1997a), (Duraes et al., 
1997b) and (Duraes et al, 1998b), (Duraes et al, 1998a), and shown here in 
Fig. 6. The leading particle can emerge from different regions of the phase 
space (cf. Figs. 1) and distribution of its momentum fraction xl is given 
by (Duraes et al, 1997a): 

F{xl) = (1 - a) dx x^"'^) ix;y = l- xl) + 

rl 

+ ^ a W dxx^'^Hx;y = i-XL), (6) 

j = 1 '2 ^min 

where a = ai + a2 is the total fraction of single diffractive (d) events (from 
the upper and lower legs in Fig. 1, respectively, both double DD and DPE 
events are neglected here) and where 



Xm.in, — Max 



itiq (Mlp + mo))^ 



(1 - xl)s' 



(7) 



with Mlp being the mass of the LP under consideration. Notice that the 
a is essentially a new parameter here, which should be of the order of the 
ratio between the total diffractive and total inelastic cross sections (Duraes 
et al, 1997a). All other parameters leading to results in Fig. 4 are the same 
as established before ^. 

We want to stress here the fact that the fair agreement with data ob- 
served in the examples shown in Figs. 4 and 5 is possible only because the 
diffraction processes have been properly incorporated in calculating the LP 
spectra (Duraes et al., 1997a). As far as the energy flow is concerned the 
IGM works extremely well (including the pionic LP not shown here but 

^In what concerns comparison with ZEUS data see also presentation of ZEUS CoUab. 
in these proceedings. Our present results differ from Fig. 4 in (Duraes et al., 1997a) where 
the preliminary ZEUS data were used instead. The only difference between the two fits 
is that whereas in (Duraes et al, 1997a) we were assuming that 30% of the LP comes 
from diffraction, now it is only 10%. 
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Figure 3. (a) Example of comparison of our LP spectra F{xl) with data from (Barton et 
al, 1983) and (Brenner et al, 2002) (Fig. 2a from (Duraes et al, 1997a)). (6) Comparison 
between our calculation and and the new data on the leading proton spectrum measured 
at HERA by the ZEUS CoUab. (c) Fits to leading J/* spectra as given by ZEUS and HI 
groups (Aid et al, 1996), (Derrick et al., 1995) and (Breitweg et al., 1997) (cf. (Duraes 
et al, 1998b) and (Duraes et al., 1998a); the fixed value of crjy^_p ~ 9 mb has been used 

and results for three different choices of p'^''* are displayed: 0.066 
solid line and 0.016 -dotted line). 



dashed line, 0.033 



discussed in (Duraes et al, 1997a)) with essentially two parameters only: 
the nonperturbative gluon-gluon cross section and the fraction of diffrac- 
tive events. At the same time, assuming the Vector Dominance Model and 
replacing impinging photon by its hadronic component (as in Fig. 3b), we 
are able to describe also the leading proton spectra observed m e — p re- 
actions. Also here the inclusion of diffractive component turns out to be 
crucial factor to get agreement with data. We can also describe fairly well 
pionic LP (not shown here, cf. (Duraes et al, 1997a)) and the observed 
differences turns out to be due to their different gluonic distributions. The 
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crucial role played by the parameter p representing the energy-momentum 
fraction of a given hadron allocated to gluons is best seen in the Fig. 3c 
example showing fit to data for leading J/^ photoproduction (Aid et al., 
1996), (Derrick et al, 1995) and (Breitweg et al, 1997). It turns out that 
(after accounting for the proper kinematics in (7) and presence of diffrac- 
tion processes as discussed above) the only parameter to which results are 
really sensitive is p = p'^/^ which, as shown in Fig. 3c, has to be astonish- 
ingly small, p'^/^ = 0.033. However, closer scrutiny shows us that this is 
exactly what could be expected from the fact that charmonium is a non- 
relativistic system and almost all its mass comes from the quark masses 
leaving therefore only small fraction, 

P ^ = ^ 0-033, (8) 

for gluons (here rUc = 1.5 GeV and Mjj^, = 3.1 GeV). 
5. Double Pomeron Exchange in the IGM 

Our latest application of the IGM discussed here will be for the DPE pro- 
cesses seen as a specific energy flow (cf. Fig. 1 d) taking place from both 
colliding particles and directed into the central region. The difference be- 
tween it and the "normal" energy flow as represented by Fig. la is that 
now the gluons involved in this process must be confined to what is usu- 




Figure 4- Our fits to two types of DPE diffractive mass distribution given by (Brandt 
et al., 2002) witfi aip-ip = 1 mb (solid lines) and 0.5 mb (daslied fines). 



ally refered to as Pomeron {IP). Such process war recently measured by 
UA8 (Brandt et al., 2002) and used (using the normal Reggeon calculus 
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arguments) to deduce the IP— IP cross section, ajp-jp. It turned out that 
using this method one gets ajp-ip which apparently depends on the pro- 
duced massAfxy- This fact was tentatively interpreted as signal of glucball 
formation (Brandt et al., 2002). However, when seen from the IGM point 
of view mentioned above, where 

1 da 2MxY dx ( MIy\ ^/,^2 2\ 

a aMxY s J xy x y xs J ^ ' 

the (Brandt et al., 2002) data can be fitted (see Fig. 4) with the same set of 
parameters as used previously to describe the SD processes (Diuacs et al, 
1997a) and (Duraes et al., 1997b) and with constant value of ajp-jp = 0.5 
mb (which is new parameter here). No glueball concept is needed here. 

6. Summary and conclusions 

The picture which is emerging from the above discussion is that the energy 
flows, which are present in all multiparticle production reactions, are ap- 
parently a kind of universal phenomenon in the following sense: they are 
all sensitive mainly to the gluonic content of the colliding projectiles (i.e., 
both to the number of gluons as given by the form of the function G{x) 
and to the amount of energy-momentum of the hadron, p, they carry and 
to the gluonic cross section which defines the actual effectiveness of the glu- 
onic component in the energy transfer phenomenon). Their sensitivity to 
other aspects of the production process (except of kinematical limits pro- 
vided by the observed masses, as illustrated in Fig. 1) is only of secondary 
importance. 

To close our arguments two other applications of IGM should be, how- 
ever, tested: whether it can also describe the final, yet unchecked energy 
flow pattern as the one provided by the Double Diffraction Dissociation 
processes and whether it can be applied in such simple form also to reac- 
tions with nuclei (first attempts were already done at the very beginning 
of the history of IGM in (Fowler et al, 1989), but they were too crude to 
be convincing at present). We plan to address these questions elsewhere. 
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